Multi-particle quantum chaos in tilted optical lattices 
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We show that, in the parameter regime of state of the art experiments on Bose Einstein Conden- 
sates loaded into optical lattices, the energy spectrum of the ID Bose-Hubbard model amended by 
a static field exhibits unambiguous signatures of quantum chaos. In the dynamics, this leads to the 
irreversible decay of Bloch oscillations. 

PACS numbers: PACS: 32.80.Pj, 03.65.-w, 03.75.Nt, 71.35.Lk 



Cold or ultracold atoms loaded into optical lattices (in 
one, two, or three dimensions) allow the experimental 
realization, with hitherto unknown precision, of a large 
variety of fascinating physical phenomena. These reach 
from quantum phase transitions 0, induced by many- 
particle interactions Q and/or disorder Q, over Bloch- 
oscillations under static forcing 0, , to quantum reso- 
nances @, quantum accelerator modes % and dynami- 
cal tunneling |3] ■ Some of these phenomena are inherited 
from the (parameter dependent) properties of the ground 
state of the atoms in the extended potential, but most of 
them are dynamical in nature, and thus imply excitations 
of the system to energetically higher lying states. 



In the present contribution, we address the spectral 
properties and dynamics of Bose Einstein Condensates 
loaded into one dimensional, tilted optical lattices, which 
we model by the Bose-Hubbard Hamiltonian (with finite 
particle number). The basic motivation is given by the 
fact that most state of the art experiments do not unam- 
biguously realize the thermodynamic limit of large parti- 
cle numbers and infinitely extended potentials, and that 
this actually is one of their decisive virtues: these systems 
allow, by systematically increasing particle number, po- 
tential extension, and, possibly, dimension, to monitor 
the convergence towards the thermodynamic limit. Fur- 
thermore, since these systems are finite, they exhibit a 
discrete rather than a continuous spectrum, with a den- 
sity of states increasing with their size. Whilst increas- 
ing the density of states leads to an increase of the time 
scale over which these systems mimic the dynamics in 
the thermodynamic limit, we shall also see that their 
basic spectral properties are uncovered once the aver- 
age level density is scaled to unity, for arbitrary systems 
size. Indeed, we will demonstrate that the spectrum of 
the tilted ID Bose Hubbard problem exhibits universal 
Wigner-Dyson statistics, a hallmark of "quantum chaos" . 
Such chaotic energy level structure implies the rapid de- 
phasing of quantum wave packets, observable, e.g., in the 
irreversible decay of Bloch oscillations. Since irreversible 
wave packet dispersion, which ultimately means decoher- 
ence, is here brought about by many particle interactions, 
we have a case of - experimentally tunable - interaction 
induced decoherence. 



The Hamiltonian of our system reads 



H 




ni 



(1) 



with a] , &i, and ni the particle creation, the particle ani- 
hilation and the number operator at site I of the lattice. 
The static held F or, more precisely, the Stark energy 
(the lattice period is set to unity), apparently destroys 
the transnational symmetry of the problem. However, the 
latter can be restored by a gauge transformation which 
amounts to changing to the interaction picture with re- 
spect to the static held, leading to 
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periodic in time with the Bloch period Tb = 1/F. Im- 
posing periodic (cyclic) boundary conditions for a hnitc 
lattice size L, we obtain the quasienergy spectrum and 
the associated eigenvectors of H by diagonalizing the 
Floquet-Bloch operator 



Cf (Tb) = exp 



T B 



H(t)dt 



(3) 



where the hat over the exponential denotes time order- 
ing. Note that, in the Fock basis {|n}}, defined by 
the iV-particle bosonic wave function constructed from 
N single particle Wannier functions U(Tb) inherits 
block diagonal structure from the translational symme- 
try of J2J, with each block labeled by a suitably gen- 
eralized quasimomcntum k of the many-particle system 
Hence, eigenstates pertaining to different values 
of k do not interact, and the corresponding eigenval- 
ues may exhibit exact degeneracies under variation of F. 
Therefore, to identify signatures of quantum chaos in the 



2 




-8 

L±J 



y/ ■ , , li 1m /M 
-■''11 Jni 



\!lMM' , / '/ J I 



/ ' //I I'' I ' >' - 

Mills ill I i I , , i , 'l sSil/h 



fil 



20 



40 



60 



1/F 



100 




0.2r 



0^ 
-0.2 



(b) 



-0.1 





Re(u) 



0.1 



0.2 



FIG. 1: Evolution of the quasienergies of the Floquet-Bloch 
operator under changes of the inverse static field 1/F. Par- 
ticle number TV = 4, lattice size L = 7, fixed quasimomentum 
k = 9], particle-particle interaction strength W = 0.032, 
hopping matrix element J = 0.038. The abundance of avoided 
crossings of any size makes good quantum numbers a useless 
concept - and identifies the present problem as a clear case 
of (many particle) quantum chaos. 



quasienergy spectrum, we have to unfold the spectrum 
such that level crossings inherited by discrete symme- 
tries of the problem do not contribute, and this simply 
means that we have to restrict to sampling spectral data 
for fixed k @. Not only is this crucial for our spec- 
tral statistics, but it also bears very pragmatic advan- 
tages in terms of the required storage space in the nu- 
merical diagonalization of U(Tb), since the dimension 
J ~ MjL = (N + L- l)\/N\(L - 1)1/ L of a single k 
block of U(Tb) is significantly smaller than the total di- 
mension AT. 

With these premises, we plot in Fig. ^ the level dy- 
namics of U(T B ), for J = 0.038 and W = 0.032 (the 
experimental values, in units of photon recoil energies, 
for rubidium atoms in optical lattices with a potential 
well depth of approx. ten photon recoils, as a function of 
the inverse static field strength. Clearly, the level struc- 
ture is "chaotic" , i.e. exhibits numerous avoided cross- 
ings on various scales (visible or unvisible on the scale 
of the figure), and it is rather obvious that good quan- 
tum numbers are no helpful concept any more (this is 
presicely the meaning of quantum chaos). Only in the 
strong field limit do the levels coalesce, since here the 
Stark energy largely dominates the hopping term and 
the particle-particle interaction in (|TJ 10]. 

Fig. |2t a) shows the cumulative level spacing distribu- 
tion 
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FIG. 2: (a) Cumulative distribution 7(s) of the quasienergy 
level spacing of the Floquet-Bloch operator (full curve), for 
F = 0.01, W = 0.032, J = 0.038, TV = 7, L = 9, and 
fixed quasimomentum k. (The dashed and dashed-dotted 
lines correspond to cumulative distributions obtained from 
Poissonian and Circular Orthogonal (COE) random matrix 
ensembles, respectively.) The numerical data unambiguously 
exhibit COE statistics, a hallmark of quantum chaos |11| . (b) 
Distribution of the real parts of the matrix elements u of the 
Floquet-Bloch operator, for the same parameter valus as in 
(a). The perfect fit by a Gaussian distribution underpins our 
identification of U(Tb) with a random unitary matrix. 



same values of J and W, though for increased system size 
N = 7, L = 9, such as to improve the spectral statistics. 
The numerical data closely follow the Wigner-Dyson dis- 
tribution for the Circular Orthogonal Ensemble (which is 
the appropriate random matrix ensemble for comparison 
with the unitary Floquet-Bloch operator [111 ]) 
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extracted from a spectrum alike the one in Fig.^ for the 



and clearly exclude Poissonian statistics (associated with 
regular spectral structure). Consistently, the real (and 
also the imaginary) parts of the matrix elements u = 
(m'|J7(TB)|n') of the Floquet-Bloch operator are Gaus- 
sian distributed, as demonstrated by Fig.|2b), what un- 
derpins our identification of U(Tb) with a random uni- 
tary matrix. We have checked that these results do nei- 
ther depend on the value of quasimomentum, nor on the 
system size (for N,L< 11, 0.5 < N/L < 1.5). Only the 
case N — L has to be excluded, due to some additional 
symmetry property which requires a separate treatment 

Let us finally address the observable consequences of 
the chaotic level structure of our many-particle Hamilto- 
nian. The average atomic momentum p(t) — tr [pp(t)] is a 
quantity readily accessible in state of the art experiments 
[4.ll2|. and determined by the single particle density ma- 
trix p(t), which can be easily deduced from the total wave 
function |\&(t)) generated by the action of U on the initial 
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FIG. 3: Irreversible decay of the Bloch oscillations of the 
average atomic momentum p(t) (saled as p — > p/NJ) for a 
weak static field F — 0.05, a filling factor n = 1 and lattice 
sizes L = 7 (top) and L = 9 (bottom). (The respective total 
Hilbert space dimensions are M = 1716 and M = 24310.) 
Comparison of both plots shows rapid convergence towards 
the thermodynamic limit N,L — > oo, N/L — const, with a 
well defined characteristic decay time r. 



condition |\&(0)). In the absence of particle-particle in- 
teractions, p(t) exhibits the well-known Bloch oscillations 
P| with Bloch frequency lob — 2nF. For non- vanishing 
W, however, the off-diagonal elements of p(t) decay on a 
characteristic time scale r, leading to irreversible decay 
of the Bloch oscillations, as born out in Fig. [3]for differ- 
ent system sizes, at fixed filling factor N/L — 1 of the 
lattice. Note that the time dependence of p(t) is virtu- 
ally identical for both choices of N and L in Fig. [JJ for 
times shorter than approx. three Bloch periods. Only 
for longer times do residual fluctuations persist, which, 
however, are quickly damped out with increasing Hilbert 
space dimension. Hence, at least as far as this specific 
dynamical observable is concerned, the thermodynamic 
limit is rapidly approached, simply due to the rapid in- 
crease of Hilbert space - and, hence, of the level density 
- with N and L. Whilst the decay of the wave packet 
is clearly due to the chaotic level structure of Fig. ^ and 
is brought about by the rapid dephasing of the various 
spectral components which constitute the initial |\l/(0)), 
we have so far no quantitative understanding of the de- 
pendence of r on W, J, and F. Clearly, this is a challeng- 
ing and equally exiting problem, since experimentalists 
will be able to enter the regime of multiparticle quantum 
chaos in a controlled way, and hence to test the genesis 
of interaction induced decoherence in the strictly unitary 
dynamics ruled by the tilted Bose-Hubbard Hamiltonian. 
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